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Many real networks are complex and have power-law vertex degree distribution, 
short diameter, and high clustering. We analyze the network model based on thresh- 
olding of the summed vertex weights, which belongs to the class of networks pro- 
posed by Caldarelli et al. [Phys. Rev. Lett. 89, 258702 (2002)]. Power-law degree 
distributions, particularly with the dynamically stable scaling exponent 2, realistic 
clustering, and short path lengths are produced for many types of weight distribu- 
tions. Thresholding mechanisms can underlie a family of real complex networks that 
is characterized by cooperativeness and the baseline scaling exponent 2. It contrasts 
with the class of growth models with preferential attachment, which is marked by 
competitiveness and baseline scaling exponent 3. 

PACS numbers: 89.75.Hc, 89.75.Da, 89.75.Fb 



I. INTRODUCTION 

Complex networks have drawn increasing interests in various disciplines. Recent studies 
have revealed that networks in the real world are far from fairly regular or totally random. 
Particularly, real networks have small average shortest path length and high clustering at 



the same time, whereas conventional graphs such as lattices, trees, or the original random 
graphs are not equipped with these properties at the same time DD 

The average path length is denoted by L, and it more or less characterizes the diameter 
of the graph. The average of the shortest path length over all the pairs of vertices defines 
L. The clustering property can be locally evaluated by the vertex-wise clustering coefficient, 
which is the number of connected triangles containing the vertex in question, normalized 
by the maximal number of possible triangles. If the vertex degree, or the number of edges 
emanating from a vertex, is k, the normalization constant becomes k(k — l)/2. The clus- 
tering coefficient C of the whole graph is the local clustering coefficient averaged over all 
the vertices. Watts and Strogatz proposed the small-world networks that simultaneously 
realize large C and small L However, the small- world networks are short of the scal- 
ing property of vertex degree distribution denoted by p{k). Indeed, not all but many real 
networks satisfy p(k) oc A; -7 typically with 2 < 7 < 3 [2J. Then, Barabasi and Albert 
(BA) developed the network model, which dynamically generates scale-free networks with 
7 = 3 0, [J. The fundamental devices in the BA model are the network growth and the 
preferential attachment; vertices are added one after another to the network, and edges are 
more prone to be connected to vertices with larger k. Various scale-free networks including 
extensions of the BA model, such as networks with dynamic edge rewiring , those with 
nonlinear preferential attachment lol, those with weights on edge s |V I , the fitness model 
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and the hierarchically and deterministically growing models (3, HO, HJJ, H2J, L 
proposed. These models largely yield more flexible values of 7, which is restricted to 3 in 
the original BA construction. Furthermore, modifications to reinforce the clustering prop- 



have been 



erty, which the BA model lacks, have also been done. A simple solution is to embed a 
triangle-generating protocol into the BA model 3] • Introduction of a node deactivation 

y designed 

0, B Q- In short, large C, small 



procedure also enhances clustering 



versions of the hierarchical models mentioned above 
L, and scale-free p(k) can be simultaneously realized by the modified BA models or by the 
hierarchical construction. Both models rely on the combined effects of network growth and 
the preferential attachment, although preferential attachment is not explicitly implemented 
in the hierarchical networks. 

Nevertheless, every network is not apparently growing. Networks can experience struc- 
tural changes that are relatively much faster than network growth or aging processes. In 



economical networks of companies, friendship networks, peer-to-peer (P2P) networks, and 
networks of computer programs, for example, it seems natural that agents change their 
connectivity without significant entries or leaves of members. Therefore, there has been a 
need for developing a nongrowing algorithm to generate realistic networks. In this regard, 
Caldarelli et al. proposed a class of networks whose connections are determined by interac- 



tions of vertices that are endowed with intrinsic weights 



considered as a type of fitness 
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The vertex weight is 



which describes the propensity of vertices 



to gain edges. It can be interpreted as money, social skills, or personal influence in social 
networks, activeness, the value of information attached to a vertex, concentration or mass of 
some ingredients in chemical or biological networks, or the vertex degree itself. Surprisingly, 
scale-free topology spontaneously emerges even with weight distributions without power laws 



118, 



20]. 



In this paper, we analyze a subclass of their model that is based on a deterministic 
thresholding mechanism. The connectivity between a pair of vertices is determined by 
whether the sum of the weights of the pair exceeds a given threshold. Actually, this class of 
networks is equivalent to the threshold graph in the graph theoretical context 2^1 , and we 
also discuss its consequence. 

Despite the stochasticity and certain continuity in the real world, thresholding that is 
more or less "hard" is often observed. Although the correspondence to our framework is not 
perfect, a common form of thresholding is that an agent on a vertex determines its action or 
state based on the number of neighbors taking a specific state. For example, propagation of 
riots, fashions, and innovations are considered to be equipped with thresholding mechanisms 
(25 1 . These phenomena have been simulated by dynamic models, such as the threshold model 



25 



, the minority games 26] , the threshold voter models, and the threshold 



for social decision 
contact processes 

In this paper, we show that a baseline power law p(k) oc /c~ 7 with 7 = 2 rather than one 
with 7 = 33] dominates this class of models and explore its cause and consequence. In 
Sec. |TIJ we follow Refs. jisl . I^J to explain the network model, and calculate fundamental 
quantities such as p(k), C, and the measure for degree correlation. The results in Sec. |H] 
are applied to various weight distribution functions in Sec. IIII[ extending the results for the 



18, 



20|. Consequently, we find that the power law is 



exponential distributions in Refs. 
observed for a wide class of weight distributions. In Sec. IIV| we argue that the power law 



with 7 = 2 is rather ubiquitous in the sense that it is a unique stable degree distribution 
when a network evolves without growth. 



II. MODEL 



Let us start with a set of n vertices V = {vi, t> 2 , . . . , v n }. As introduced in Refs. |l8l. l20j. 
we assign to each Vi (1 < i < n) a weight W{ G R that is taken randomly and independently 



distributed as specified by a given probability distribution function 
quantifies the potential for the vertex to be linked to other vertices 
the weight permits additive operation. Actually, the multipliable weights 
can be easily reduced to the additive ones by taking the logarithm of w. Let 



w) on R. The weight 
la. We assume that 
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F(w) 



f{w')dw' 



(1) 



and lim^^oo F(w) 



be the cumulative distribution function, satisfying lim„,_ > _ 00 F(w) 

1. The set of edges E is defined by the thresholding rule with threshold 6: E = { (_ Vj , vj ) ; Wi + 



Wj > 6,i 7^ j}. We focus on this specific case of more general framework [la, |2Qj. This 
renders more mathematical analysis and allows us to explore the consequence of vertex 
interactions based on intrinsic weights. The degree distribution p(k), where < k < n is the 
vertex degree, is readily calculated with the use of continuum approximation corresponding 
to the thermodynamic limit (n — ► oo). However, we confine ourselves to a finite n, and the 
limit n — > oo should be understood as approximation. Putting the upper limit of k equal to 
n instead of n — 1, we obtain 



k = n f(w')dw' = n [1 - F (6 - w)} , (0 < k < n) 
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and 



p(k) = /O)^ 



1 



0) 



(2) 



(3) 
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Because of the one-to-one correspondence between k and w represented by Eq. (j2J), 
the vertex-wise cluster coefficient depends only on k, which simplifies the analysis. We 
denote it by C(k\ and the scaling law C(k) oc k~ x is often observed in real and mod- 



eled networks [1( 
is given by C = 



12 



13, 121 



The clustering coefficient of the entire graph 
S-cm^ To calculate C(k) let us consider a vertex v 



with degree k = n[l — F(0 — w)\. The density of the number of neighbors with degree 
k' — n [1 — F(0 — w')] becomes f(w') if w' > — w and otherwise. With such a neighbor 
denoted by v', the number of connected triangles comprising v, v' and another neighbor of 
v is obtained as follows. When k' > k, a new neighbor of v is also a neighbor of v' because 
w' > w. The number of triangles in this case is 

(n - 2) [1 - F(0 -w)]^n[l- F(6 - w)] = k. (4) 

When k' < k, we have 



nf(w")dw" = n[l- F(6 - w')] (5) 

triangles. We obtain C{k) by the sum of Eqs. (JU) and (0) that is weighted by the degree dis- 
tribution. The normalization is given by dividing it by k(k — l)/2 and by another factor of 2, 
as each triangle is counted twice. Consequently, we have for w > 8/2, or k > n [1 — F (0/2)] 

°« = ^M^{^ fc/(w ' Mw ' + /.l" [l - F(s - TO ' )l/(w ' ) ' <, "'} 

-/* fl--)p(A;')^)/(A:/n) 2 . (6) 

When w < ^/2 or < n [1 — F (^/2)], we simply end up with 

C7(fc) = I. (7) 

The vertices with C(k) = 1 forms the peripheral part of the network that is connected to the 
cliquish core with smaller C(k), as schematically depicted in Fig. ^ The core consists of the 
vertices with w > 0/2, and it is similar to the winner-take- all phenomenon found in growing 
network models"^ 0]. However, more than a single winners are allowed in this model. This 
core-peripheral separability is a rigorous property of the threshold graph 2^|. It is also 
consistent with the property of real networks that C(k) saturates for small k |3, [3] and 



that vertices in the core are densely connected 

Regarding the network size, real networks are small with L proportional to Inn or even 
less QjsHsj]. Our network has L < 2 as far as it is free from isolated vertices. This is because 
any pair of vertices can be connected by a path of length 2 passing a vertex with a sufficiently 
large weight is in the cliquish part (see Fig. [l]). One might ascribe this ultrasmallness to 



the fact that the mean degree is of the order of n, indicating too many edges per vertex. 
However, the mean degree can be kept finite by scaling 9 according to the increase in n, as 
discussed in Sec. IIII Al It turns out that this modification does not change L. 

The correlation between the degrees of adjacent vertices also characterizes networks 
Actually, degree correlation can be positive or negative depending on the type of network, 
as measured for real data using the degree of assortativity |30|. Here we explain a simpler 
quantity to gain insight into the degree correlation P, [lg, Il7| . which was first analyzed in 
Ref. Q for the present type of network model. We denote by k 2 the sum of the degrees of 
v 's neighbors, which has degree k. If the degree is uncorrelated, kijk, or the average degree 
of the neighbors, is independent of v or k. For the threshold model, we derive 



roo 

k 2 = / nf(w')n [1 - F (0 - w')) dw' 
Je-w 



n 2 



--f , (l--)p(k')dk'\- (8) 



Accordingly, 
which generally depends on k. 



h = n J i _ I ( n _ k')p(k')dk'\ , (9) 

k I k Jn[l-F(6-F-i(l-k/n))] J 



III. EXAMPLES 



In this section, we calculate the quantities introduced in Sec. |H] for some weight distri- 
butions f(w). 

A. Exponential distribution 

Let us begin with recapitulating the case of exponential weight distribution with deter- 
ministic thresholding [3,0]. We set 

f(w) = Xe- Xw (0 < w). (10) 

Assuming 6 > so that the generated networks are not trivial, for w < 8, combining Eqs. (jHJ) 
and (fTUj) results in 

p(Jfe) = [ne- xe <k<n). (11) 

k 



A total of nf(w)dw = 
In general, this sort 



The scale-free distribution p(k) oc k~ 2 appears from random weights whose distribution has 
nothing to do with power law, which is a main claim of Ref. Hi 
ne~ xe vertices are condensated at k = n and form a core 
of condensation occurs when f(w) has a lower cutoff. However, the number of vertices 
with k = n can be made arbitrarily small by setting large X9, and this feature is not so 
essential. We numerically simulate a network with n = 50000, which is fixed throughout 
the paper, A = 1, and 9 = 10. Setting A = 1 does not cause the loss of generality because 
only the multiple of A and 9 appears in Eq. (|TT|) and the following quantities [Eqs. (fT2"j) and 
« Figure HaUhows that _a> results (crosses) are pred.cted by Eq. CO (lines) 



sufficiently well jj, 3- 
C{k) = < 



In regard to clustering, Eq. © yields 



(ne 



-xe 



<k< ne- xe / 2 ) 



fc 2e 



-xe 



1 + X9 + 2 In |) {ne- xe l 2 < k < n) 



(12) 



which agrees with the numerical results in Fig.|2fb) (crosses) (originally derived in Ref. 
Equation fT2|) shows that C{k) nearly decays according to the power law with exponent 
2. However, analysis of real networks, such as metabolic networks Q], actor networks, 
semantic networks, world wide webs, and the Internet Jl3| , suggests C(k) oc k~ l , which is 



also supported by some models 



|28| . For this particular example, we 



have a larger scaling exponent. Actually the general power-law form C(k) oc k 7 ' (7' > 0) is 
also reported in model studies and in data analysis p. The clustering coefficient 

of the whole network is 



C 



C(k)p(k)dk 



-xe/2 



9 



5 + 3X9 



-2xe 



(13) 



In real networks, the average vertex degree denoted by (k) is independent of n on a large 
scale 
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. Since 



(k) 



-xe 



n + X9) 



(14) 



finite (k) is maintained by setting 9 = A -1 Inn. With this scaling of 9, our model produces 
a finite value of C that does not decay to in the limit n — > 00. This result agrees with real 
data 2:], and C is actually nonvanishing for more general f{w). 

Using Eq. 0, the average degree of neighbors becomes, as shown in Ref. j^ . 

k\ 



k 



n 2 e -xe 



k 



1 + A# + In • 



n 



(15) 



Since Eq. ()15j) is decreasing in k, the network is disassortative |30J] with negative degree 



correlation, which is a property shared bv some scale- free network models |5|, |16j, ^| and 
some social and biological networks 0, 0, 0|. More intuitively, disassortativity of our 
model is a natural consequence of the core-peripheral structure. Equation (fTK|) and Fig.Efc) 
(crosses and solid lines) show that there is an approximate scaling relation k 2 /k oc k^ 1 [17]. 



B. Logistic distribution 

The logistic distribution is often used, for example, in statistics and economics. Except 
the discrepancy in the asymptotic behavior, it can be used as a more or less appropriate 
substitute for the Gaussian distribution. The logistic distribution is more tractable because 
it has an analytic form of F(w). For a given f3 > 0, the logistic distribution is defined by 

f( w ) = —Bl (16) 

M ; (1 + e-^) 2 V ; 



With 



F-^) = -im(i-i), (is) 

k = i + e 0(w-o) ' ( 19 ) 

w = e + \\n— k —. (20) 
p n — k 



applied to Eq. (jSJ), we obtain 



p{k) = —-^ — 2 (0<k<n). (21) 

k 2 (l + e-^2=*) 

The power law p(k) oc k~ 2 is again manifested as k approaches n. If 1 [(n — k)/k]e~@ e , k 
is relatively small, and p(k) = ne@ /{n — k) 2 does not strongly depend on k. The crossover 



from this regime to the power-law regime, which is found in real data 12, [l3( and derived 
by scaling ansatz theory Q], occurs around 1 = \{n — k)jk\e~^ b \ or k = [n/(e^ e + 1)]. A 
larger value of (38 provides a wider range of k in which the power law holds. In this range, 
Eq. (fi7)|) is of course approximated by the exponential distribution represented by Eq. (fTUj) 
with (3 = A. For k small relative to n/i^e 139 + 1), f(w), with a correspondingly small w, does 



not decay exponentially or even monotonically. Therefore, when k is small, the number of 
vertices with degree k is not large enough to support the power law. 

We compare in Figs. Efa) the numerical results for 9 = 6 (open squares) and 9 = 10 
(open circles) with the corresponding theoretical results in Eq. (12 lj) (dotted lines). We have 
set (3 = 1 without losing generality for the same reason as in Sec. IIII Al The effect of 9 on 
the position of crossover is clear in the figure. Since the integrals in Eqs. © and Q cannot 
be explicitly calculated, numerically evaluated C(k) and k 2 /k with the same parameter 
values are shown in Figs. Efb) and|2fc), respectively. Similar to the case of the exponential 
distribution, the crossover from the plateau to the power law is observed for C(k) with the 
same scaling exponent C(k) oc k~ 2 . Also with regard to the degree correlation, k 2 /k oc k~ l 
approximately holds except for small k. 



The Gaussian distribution can be a standard null hypothesis on the weight distribution. 
Since it does not have the analytical form of F^ 1 {x), we perform straightforward numerical 
simulations with 9 = 6 and 9 = 10 to examine p(k), C(k), and k 2 /k. The Gaussian 
distribution is assumed to have mean and standard deviation 1.7 to roughly approximate 
the logistic distribution with (3=1, which has been used in Sec. IIII Bl In spite of different 
asymptotic decay rates of f(w), Fig. |2] indicates that p(k), C{k), and k 2 /k for the Gaussian 
distribution do not differ so much from those for the logistic distribution, disregarding the 
crossover points. This implies a rather universal existence of power law behavior, which is 
discussed in more detail in Sec. IIVI 



The Pareto distribution, which is equipped with an inherent power law, is often observed 
in, for example, distributions of capitals and company sizes jsj]]. It is defined by 



where a > and Wq > 0. Nontrivial networks form if we choose 9 > 2u> . We obtain 



C. 



Gaussian distribution 



D. Pareto distribution 




(22) 




(23) 



When w < 9 — wq, it is straightforward to derive 

w 



k = n 
w = 9 
p(k) 



9-w 

n \ l/a 

k) 

n l/a 



n 



w 



9 -wq 



< k < n 



Wo, 



(J_ k l/a 



n 



l/a 



a+l 



A total of 



nj[w)aw = n 

8-w \9 — Wo 



(25) 
(26) 
(27) 

(28) 



vertices with w > 6 — wo are condensated at k — n. When n{wo/6) a <C k < n, p(k) can be 
approximated by a power law 

a+l 



p(k) 



Wo 

e 



n l/a£-(a+l)/a_ 



(29) 



By modulating a, we can produce a scale-free p(k) with arbitrary 7 = (a + l)/a > 1. An 
observed 7 in turn serves to estimate a and f(w), which may underly, for example, fractal 
dynamics of economical quantities, as well as network formation. The scaling exponent 
for p(k) differs from that for f(w), and a faster decay of f(w) with a larger a yields a 
slower decay of p(k). Numerical results for p(k) are shown in Fig. Efa) with wo = 1. We 
set (a, 9) = (0.5,100) (open squares) and (0.5,500) (open circles), yielding 7 = 3, and 
(a, 9) = (1, 100) (closed squares) and (1, 500) (closed circles), resulting in 7 = 2. The results 
are consistent with the theoretical prediction based on Eq. (|2Tj) (solid lines) and also with 
Eq. (JH (dotted lines). 

and (|2*7jl into Eqs. ©, ©? and Q, respectively, yields 

„2a-l 



Substituting Eqs. 
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k > n 



V 9 



C{k) 



n 



Wq 

9 -Wq 



< k < n 



(2w, 



(30) 
(31) 



and 



k 1= [_n_a (woy ft** X ^zM^l dx 
k I k V 6 J Je/[e-w (n/k) 1 / a ] (x — l) aH 1 



2+1 

where we set y = x -1 . The integral in Eq. ({HQ]) is nonnegative and does not depend so much 
on k when k tends large. Therefore, C(k) oc k~ l is expected based on the first term, which 
is consistent with the numerical results for k larger than the crossover value k = n (2wo/9) a 
[Fig- Eft 1 )]. The scaling law C(k) oc A; -1 , as opposed to C(k) oc k~ 2 for the exponential f(w), 
rather agrees with real data jl3 fl^l] . 

Similarly, Eq. (jH2j) and the simulation results shown in Fig. |^c) suggest k^jk oc k~ x for a 
sufficiently large k. Equations (fHUj) and (jH^jl show that the scaling exponents of both C(k) 
and /C2//C do not depend on 7 or a. 



E. Cauchy distribution 

For the Cauchy distribution 

f(nn\ = 

7r(l + W 2 



f(w) = n ] 9 . (weR), (33) 



we obtain 



F-\x) = tan|(2x-l), (34) 
71 ( 2k\ 

w = - tan - ( 1 J , (35) 



1 + tan 2 § (l 



2A- 



p(*) = \ n ' 2 (0<k< n). (36) 

»l+(0-tanf (l-f)) 

Numerically obtained p(fc), C(k), and ^2/^ together with Eq. ()36|) are shown in Fig. 0] 
for 6 1 = 100 (open squares) and for 9 = 500 (open circles). According to Eq. (jSHJ), the 
monotonicity of p(k) is marred because p(0) = p(n) = 1/n. A particular choice of 9 = 
even gives rise to the uniform p(k). For general 9, however, p(k) has the unique maximum 
and minimum between k = and k = n as shown in Fig. IHa) by solid lines. Existence of the 
characteristic vertex degree corresponding to the peak of p(k) is a feature shared by random, 
regular and small-world networks 

HQ. 

The peak appears because of the unimodality of 



f{w), which yields the plateaus of p(k) in the case of the logistic and Gaussian distributions 
(see Sees. IIII Bl and IIII Cj) . Nonetheless, as for the Pareto distribution with a = 1, which 
has the same asymptotics f(w) oc w~ 2 as the Cauchy distribution, approximate power laws 
with 7 = 2 are observed for intermediate values of k. 

The one-sided Cauchy distribution on the half line facilitates a fairer comparison with 
the Pareto distribution that also has the lower cutoff of w. We define the one-sided Cauchy 
distribution by 

fM= 2 (w>0). (37) 

7r(l + W Z ) 

Then it holds that F~ 1 (x) = tan(7r/2)x and 



2 l + tan 2 ^(l-M 
p(k) = 2 \ nJ (0<k<n). (38) 

Figure EJa) shows the numerical results for p(k) with 9 = 100 (closed squares) and 9 = 500 
(closed circles), accompanied by the prediction by Eq. (jHHjl (dotted lines). The approximate 
power law holds even for k close to n, which contrasts with the case of the standard Cauchy 
distribution. The behavior of C(k) and fo/fc shown in Figs. Efb) and|3Jc), respectively, 
resembles that for the Pareto and standard Cauchy distributions. In Sees. IIII Al and IIII Bl 
we have inspected consequences of using exponential f(w) and logistic f(w). Including the 
comparison between the Pareto and Cauchy distributions examined here, effect of a lower 
cutoff of f{w) does not seem so prominent. 

IV. WHY POWER LAW WITH 7 = 2? 

The power law of p(k) with 7 = 2 seems universal for thresholding mechanisms not only 
because a wide class of f{w) generates it but also owing to its stability. To be more specific, 
weights can be the vertex degrees themselves, as is implied by the BA model. Indeed, k 



22 1 . Then let us iterate our construction 



can represent how central or influential a node is 
algorithm to simulate an evolving but not growing network with dynamic f{w) and p(k). 
Let us simulate a dynamical network with n = 50000. Initially, w is uniformly distributed 
on [0, 1], and the thresholding algorithm determines k. Then we set w = k/n + £, where £ is 
the Gaussian white noise with standard deviation a = 0.2, and iterate the dynamics. The 
numerical results are shown in Fig. [S]with 9 = 1. In the early stages [crosses in Fig. Efa)], 



k is distributed more or less uniformly since a uniform f(w) yields a uniform p(k) just 
accompanied by a singularity at k = or k = n, which is easily checked with Eq. (J3J). As the 
iteration goes on, however, p(k) converges to a power law with 7 = 2. Similarly, C(k) oc /c _1 
and k2/k oc k" 1 are eventually realized as shown in Figs. Efb) and 03(c). Although 9 too 
far from 1 or excessively small a results in a complete or totally disconnected network, 
p(k) oc k~ 2 emerges robustly against changes in 9 > 1 and a > 0.15 unless noise is not 
extremely large. 

By the analogy of the Pareto case, if a power law with 7 7^ 2 is obtained, 7 will be 
transformed by the map 7 = 0+1—^7 = (a + l)/a, namely, 7 — > 7/(7 — 1). This map 
has a unique positive fixed point 7 = 2. Actually, the map is neutrally stable at 7 = 2 with 
eigenvalue —1, which implies oscillation. This argument does not directly support but may 
underlie the emergence of p(k) oc k~ 2 . In addition, p(k) converges to p(k) oc k~ 2 stably with 
respect to the choice of initial distribution f(w). It is in a striking contrast with the case 
of competitive growing networks with vertex weights, which generate p(k) oc k -1 only for 
a Hmited e.ass of weight dotations Q. More broadly, general cooperative networks in 
which interactions of multiple vertices leads to interconnection [18] may have stable power 
laws, possibly with 7/2. As an example, we can assume w of the next generation to be 
proportional to k x {x > 0), which is often the case in real networks |22j. In this case, it is 
easy to show that the neutrally stable fixed point of the abovementioned map is 7 = 1 + y/x. 
Our model, which yields 7 = 2, sets a baseline example of this class. 

In terms of clustering, we have observed the discrepancy in the scaling law C(k) oc k~ 2 
for the exponential, logistic, and Gaussian distributions with C(k) oc k^ 1 for the Pareto, 
Cauchy, and one-sided Cauchy distributions. The Pareto distribution with a = 1 results 
in p(k) oc k~ 2 , which coincides with the scaling law for the exponential type of f(w). It 
means that the consistency in p(k) for different choices of f(w) does not necessarily mean 
the consistency in network structure. The difference between the exponential tail and the 
power-law tail of f(w) is likely to cause qualitative discrepancy in C{k). On the other hand, 
the degree correlation behaves similarly in all the examined cases, namely, kijk oc k^ 1 . 
From a dynamical point of view, C(k) evolves from a general form to C(k) oc k~ 1 , which 



is realistic 



131 ] . This scenario matches the simulated dynamics shown in Fig. 03(b). In 



this case, C(k) is initially just large for most vertices and converges to C(k) oc k 1 , which 
reflects the eventual separation of the network into the core and the peripheral part. 



Boosted by the original BA models, the power law of p(k) with 7 = 3 has been pronounced 
in the first p.ace QQ. Moreover, in perception and contact processes on scale-free networks, 
the critical value of the infection rate is extinguished if 7 < 3 32J. These results suggest 
relevance of the power law with 7 = 3. However, real scale-free networks have more dispersed 



values of - 
2 and 00 I 



2|, and many models have been 



a33 



10 



11 



12 



13 



vertex weights (also see Refs. 



18 



proposed so that 7 is tunable somewhere between 
. In contrast, the present model with intrinsic 
20j l) and another type of thresholding model 3^| broadly 



14. 



16 



22 



yield 7 = 2. We speculate that 7 = 2 is another general law. In the parameter space of 7, 
2 as well as 7 = 3 often emerges as phase transition points of network characteristics 
2Ji|22]. The 7 = 2 law may be common to cooperative models such as those with 



thresholding, while the 7 = 3 law underlies the competitive models represented by network 
growth with preferential attachment. Actually, many real networks in the scale- free regime 

Q 

have 7 close to 2 rather than to 3 [2j. Some networks such as the world wide web, e-mail 



networks, language networks, and ecological networks have 7 even smaller than 2 |24|. Some 
of these observations can be understood as small deviations from our 7 = 2 law, which 
■nay be explained by pr 0£ ,r .edification of the m ode. HQ, for e— by introdncing 



stochastic thresholding 
interactions. 



21 



, nonlinear relations between k and w |22j, or many-body 



V. CONCLUSIONS 

We have shown that the thresholding model, which is in the class of networks with 
intrinsic vertex weights, generates scale-free networks with 7 = 2, large C, and small L for a 
broad choice of weight distributions. Even if we start with an arbitrary weight distribution, 
p(k) oc k~ 2 and C(k) cx k~ x are finally obtained. The competitive mechanisms, such as 
network growth with preferential attachment or hierarchical structure, are not mandatory 
for generating realistic networks ^J. The cooperative thresholding mechanisms also result 
in desired properties rather generally, and they yield somewhat different characteristics from 
those of growing types of networks. In addition, they allow plausible physical interpretations, 
have core-peripheral structure, are equipped with inhomogeneity as in real networks, and 
facilitate analytical calculations 
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Figure captions 



Figure 1: Schematic diagram of the threshold model. 

Figure 2: Numerical results for (a) p(k), (b) C(k), and (c) k 2 /k using the networks of 
size n = 50000 generated by thresholding. The weight functions are taken to be exponential 
with A = 1, 9 — 10 (crosses), logistic with (3 = 1, 9 = 6 (open squares) and /3 = 1, 9 — 10 
(open circles), Gaussian (mean and standard deviation 1.7) with 9 = 6 (closed squares), 
and Gaussian with 9 = 10 (closed circles). The theoretical predictions are shown for the 
exponential distribution (solid lines) and the logistic distributions (dotted lines). 

Figure 3: Numerical results for (a) p(k), (b) C(k), and (c) k 2 /k for the Pareto weight 
distributions with n = 50000 and w = 1. We set a = 0.5, 9 = 100 (open squares), a = 0.5, 
9 = 500 (open circles), a = 1, 9 = 100 (closed squares), and a = 1, 9 = 500 (closed circles). 
In (a), p(k) estimated by Eq. (|2*7j) and the by power law approximation in Eq. (|29jl are also 
shown with solid lines and dotted lines, respectively. 

Figure 4: Numerical results for (a) p(k), (b) C(k), and (c) k 2 /k for the Cauchy weight 
distribution with n = 50000 with 9 = 100 (open squares), 9 = 500 (open circles), and the 
one-sided Cauchy distribution with 9 = 100 (closed squares), and 9 = 500 (closed circles). 
In (a), the analytically estimated p(k) for the Cauchy distribution [Eq. (|36))] and for the one- 
sided Cauchy distribution [Eq. ([38)1] are also shown by solid and dotted lines, respectively. 

Figure 5: The evolution of (a) p(k), (b) C(k), and (c) k 2 /k with the repetitive thresh- 
olding. We set n = 50000 and 9=1. The data shown are those after 1 (crosses), 8 (open 
squares), 10 (closed squares), 12 (circles), and 15 (triangles) rounds. 
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